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The Pb/Si(lll) thin films were simulated within the density functional theory (DFT). The well- 
known Perdew-Burke-Ernzerhof (PBE) version of the generalized gradient approximation (GGA) 
and its recent nonempirical successor Wu-Cohen (WC) issue were used to estimate the exchange- 
correlation functional. Lattice parameters were calculated for Bulk of the Pb and Si compounds to 
obtain more reliable lattice mismatch at the interface to be consistent with our used full-potential 
method of calculations. The WC-GGA result predicts the lattice constants of the Pb and Si com- 
pounds better than the GGA when compared with experiment. We have found that the spin- 
orbit coupling (SOC) does not significantly influence the results. Our finding is in agreement with 
the recent observation of the Rashba-type spin-orbit splitting of quantum well states in ultrathin 
Pb/Si(f f f ) films. Our result shows, in agreement with experiment, that the top site (Tl) is the most 
stable phase. A combination of tight o and feeble 7r bonds has been found at the interface, which 
results in a covalent Pb-Si bond. Our calculated electric field gradient (EFG) predicts quantum 
size effects (QSE) with respect to the number of deposited Pb layers on the Si substrate. The QSE 
prediction shows that the EFG dramatically drops on going from first to second layer. The EFG 
calculation shows that this system is not an ideal paradigm to freestanding films. 

PACS numbers: 71.15.Mb, 73.20.At, 76.80.+y 



I. INTRODUCTION 



The Pb/Si(lll) thin films can be considered as an ideal 
interface, since it is an unreactivc metal-semiconductor 
system which forms perfectly an abrupt interface with- 
out intermixing. [l[ Intensive experimental 0, 0| and 
theoretical [1, [1] works performed on this system show 
that it is a subject of intense technological and scientific 
interest with a wide range of applications. Several struc- 
tures were proposed experimentally for this system, [l], Q 
The structures vary with their circumstances. The varia- 
tions form a complicated phase diagram. The stability of 
the proposed structural models were studied theoretically 
employing pseudopotential calculations. Q The complex 
phase diagram shows that there are 4 structures with a 
simple (lxl)-unit cell. The experimentally observed^ 
simple unit cell provides a suitable case to perform ac- 
curate atomic relaxations within an all-electron calcula- 
tion reliably. Our full-potential surface calculations were 
then devoted to study these (lxl) structures with the 
concentration on their most stable phase. Quantum size 
effects (QSEs) in nanostructures, e.g., a two-dimensional 
(2D) ultrathin metal films, appear as dramatic oscilla- 
tions in many physical properties upon variation of film 
thickness. 0, The motion of electrons in the film plane 
of such a 2D-systcm is essentially free, whereas the elec- 
trons are confined in the normal direction to the film 
surface, which leads to the quantized electronic states, 
i.e., quantum well states (QWS).[ll[ The oscillation pe- 
riod was found to be 7-monolayer-height-islands of lead 
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grown on top of a silicon substrate for the electron in- 
tensity curve as a function of the normal component of 
the electron momentum transfer. [8[ Vazquez de Parga 
and coworkers observed the QSEs for a number of sur- 
face properties, i.e., the surface roughening temperature, 
work function, chemical reactivity, or the surface dif- 
fusion barrier, in the Pb/Si(lll) and Pb/Cu(lll) thin 
films by Scanning Tunneling Microscopy/Spectroscopy 
(STM/STS).[| The step height of layer N is related to 
the electron spillage length into the vacuum which de- 
pends on the DOS at the Fermi level. @ Quantum size 
effects (QSEs) were found for the total energies and en- 
ergy differences in freestanding Pb(lll) thin films by us- 
ing pseudopotential method as embodied in the CASTEP 
code.[l(| Influence of the Cu(lll) substrate on the geom- 
etry structure of the Pb(lll) layers was also indirectly 
taken into account by solely reducing in-plane lattice 
constant to impose 3.3% lateral compression on all the 
Pb(lll) layers to mimic the role of Cu(lll) substrate. [l(| 
The effects of Cu(lll) substrate were found to be sub- 
stantial, as the step heights were calculated to be in 
more agreement with experiment within the later in-pane 
strain, [ijj] Wei and Chou theoretically, employing ultra- 
soft pseudopotential calculations, studied the quantum 
size effects (QSEs) in the clean surface of Pb(lll).[l2T| 
Recently P. S. Kirchmann et al. experimentally ob- 
served quantum size effects (QSEs) in the Pb/Si(lll).[l]| 
They[13[ found their experimental results for the ac- 
tual Pb/Si(lll) slab to be in agreement with the pseu- 
dopotential results [l2j for the hypothetical freestanding 
Pb(lll) slab. Dil et al.[lj] comparing their Pb/C full- 
potential calculations with the Pb/Si(lll) experimental 
results of the others [l^ reported that unlike for Pb on 
graphite, the Pb overlayer lattice structure is influenced 
by that of the Si(lll). Here we then intend to, going 
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beyond the hypothetical freestanding Pb(lll) approxi- 
mation, quantitatively assess the effects of electronic and 
crystalline structures of the underlying Si substrate on 
the deposited Pb films. For this reason we inspect the 
electronic structures at the interface of the Pb/Si(lll) 
slab. From our electronic structure calculations, we show 
that there is a mixed state composed of weak 7r and 
strong a bonds between Pb and Si at the interface which 
results in a strong Pb-Si covalent bond. We aim to in- 
vestigate whether the substrate can influence the results. 
In order to accomplish the investigation, we have calcu- 
lated the work function, energy differences and surface 
formation energy as a function of number of Pb layers 
for the (1 x 1)-Pb/Si(lll). Our result shows that the 
effect of Si substrate can be of significant importance 
for some physical quantities depending on their sensi- 
tivity to the valence electron charge density. We have 
observed the quantum size effects (QSEs) for the (1 x 1)- 
Pb/Si(lll) thin films. The goal of this work more trans- 
parently is achieved by presenting a physical interpreta- 
tion for our QSE calculations for this system in the elec- 
tric field gradient (EFG) as an extremely sensitive quan- 
tity to the valence electron charge density distribution. 
We would also examine the effect of spin-orbit coupling 
(SOC) on our thin films. The result, in agreement with 
experiment. fill] shows that the SOC has a minor effect 
for the Pb/Si(lll) thin films. Low temperature \fi x \/3, 
3x3 and \fl x \[Z phases of Pb/Si(lll) were theoreti- 
cally studied employing pseudopotential method. [l?!] It 
was shown that there could be a discrepancy in predict- 
ing the ground state of Pb/Si(lll) system between lo- 
cal density approximation (LDA) and generalized gradi- 
ent approximations (GG A1. 11711 Another objective of our 
density functional theorv|la.[l9j (DFT) calculations is to 
go beyond the LDA and GGA by using the noncmpiri- 
cal GGA recently proposed by Wu-Cohen (WC) for the 
exchange-correlation functional, [20j | which is expected to 
improve metal surface formation energy and lattice con- 
stant calculations [10, [H| compared to the Perdew-Burke- 
Ernzerhof GGA (PBE-GGA).[H 



II. COMPUTATIONAL DETAILS 

State-of-the-art calculations were performed within the 
density functional the qrvfi~8l . [l9j (DFT) as implemented 
in the WIEN2k code.[23| The augmented plane waves 
plus local orbital (APW+lo) method [H, [2al has been 
used for solving the Kohn-Sham equations [19( employ- 
ing the latest version of the generalized gradient approx- 
imation (GGA), i.e., Wu-Cohen (WC), for the exchange- 
correlation functional ppj. We set the Muffin-tin radii to 
Rmt = 2.2 bohr for the Si and to Rmt = 2.5 bohr for 
the Pb atoms. The expansion of the wave functions and 
charge densities were cut off by the RuT^-max = 7.5 and 
G m ax = 14 parameters, respectively. The full relaxations 
were performed with the criterion of 1 mRy/bohr on the 
exerted forces. The relativistic effects were taken into 



TABLE I: Lattice parameter, bulk modulus and pressure 
derivative of bulk modulus of Si (Pb) together with the 
method of calculations and exchange-correlation (XC) func- 
tional. 



Method XC-Potential a(bohr) B(Gpa) 



APW+lo 

APW+lo 

APW+lo 

PPW 

PPW 

Exp. d(e) 



GGA-WC 
GGA-PBE 
LDA 
LDA a(f 

GGA c(6) 



10.28 (9.31) 
10.34 (9.49) 
10.22 (9.21) 
10.17 (9.16) 
10.32 (9.56) 
10.26 (9.36) 



93 (47.3) 
88 (42.2) 
96 (53.1) 
87 (52.6) 
87 (39.3) 
98 (43.2) 



B' 
4.23 (2.9) 
4.16 (5.0) 
4.37 (3.7) 
3.82 (5.2) 
- (4-7) 
4.02(4.9) 



a Reference 
b Reference 
c Reference 
d Reference 
e Reference 



account by including the spin-orbit coupling (SOC) in a 
second variational procedure. A set of 22 x 22 x 1 special 
k points has been used for integrations over the Brillouin 
zone of the lxl surface cell. 



III. BULK STRUCTURAL PROPERTIES 

Before going through the interface and surface calcu- 
lations in detail, it seems advisable to first shortly re- 
port on the bulk structural properties. In this case, 
the accuracy of our calculations could be first assessed, 
and second more accurate lattice parameters would be 
also obtained to simulate the substrate with an actual 
and reliable lattice mismatch at the interface. Fur- 
thermore one aims to get more insight into the var- 
ious exchange-correlation functionals, specially the re- 
cent proposed WC-GGA[2C| potential. The lattice pa- 
rameters, bulk moduli, and the pressure derivative of 
the bulk moduli were calculated for the bulks of the Si 
and Pb compounds. The calculations were performed 
within the LDA[I|, PBE-GGA[13 and WC-GGA0 
functionals by fitting the total energy as a function of 
volume with the Birch equation of state. [27j The cur- 
rent APW+lo and previous pseudopotential plane wave 
(PPW) results [H, [29|, [3(| along with the experimental 
data (3l|, [H| are presented in Tab. m Our APW+lo re- 
sult, compared to the experiment pll. HH, shows that, 
for both the Si and Pb cases, the PBE-GGA[m overesti- 
mates the lattice constants and bulk modulii, whereas the 
LDAprjj underestimates them. This is consistent with 
the general features of the PBE-GGA[i3l and LDA[H 
functionals for solid states. In particular, as shown in 
Tab. HI our values are in complete accord with the pre- 
vious PPW results. [H, M, Hfl] The WC-GGA[H pre- 
dicts the lattice constant of the Pb compound better 
than the LDA [HI and GGA regardless of the APW+lo 
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TABLE II: The work function, (/>Si (eV), and surface forma- 
tion energy, E' Si (eV /bohr 2 ), per unit area for the number of 
Si(lll) bilayers (N Si ). 





Sl (eV) 


E f si (eV/bohr 2 ) 


4 


4.8484 


1.3037 


6 


4.7179 


1.3281 


8 


4.7295 


1.3274 


10 


4.7299 


1.3270 



or PPW method, while our PBE-GGA0 remains su- 
perior for the prediction of its bulk modulus. This can 
also be the case for the lattice parameter of the Si com- 
pound. From the result presented in Tab. [I] one can cau- 
tiously state that the WC-GGAj2p| functional improves 
the lattice parameters over both the LDA[26j and PBE- 
GGA[22|. This statement is in agreement with the pre- 
diction of Wu and Cohen for the most of their studied 
compounds [20l] . Our results show that lattice param- 
eters calculated from WC-GGA[20j approximation arc 
very close to experiment |3lT IHJ. The calculated lattice 
parameter of the bulk Si, i.e., 10.28 a.u., will be used in 
the forthcoming Sec. II VI for simulating this compound as 
our substrate to mimic the effects of the lattice mismatch 
on the valence bands at the interface of the Pb/Si(lll). 
Finally, keeping in mind that the value obtained for the 
pressure derivative of the bulk modulus depends on the 
range of fitting and as a result quantitative analysis may 
not be completely reliable, our result shows that the 
PBE-GGA[22[ yields better pressure derivatives of the 
bulk moduli for the Pb and Si compounds than LDA[26| 
and WC-GGA0. 



IV. SLAB DETERMINATION AND 
OPTIMIZATION 

The silicon substrate has been simulated using the cal- 
culated WC-GGA[20| lattice parameter in the previous 
Sec. IIIII along the (111) crystallographic direction of the 
Fd3m space group. We have found that the Si(lll) sub- 
strate, as shown in Fig. [TJa), can be properly modeled 
by stacking 8 Si bilayers, i.e., 16 Si monolayers. Surface 
formation energy, E Si , per unit area and work function, 
cf>Si, were calculated for the various number of Si(lll) 
bilayers, N$i, and the results are listed in Tab. El O ur 
result shows that the changes of the work function and 
the surface formation energy are less than 0.001 eV by 
adding more bilayers. The reliability of the obtained 8 Si 
bilayers has been further investigated by comparing the 
bulk and surface electronic structures of the silicon. The 
total densities of states (DOSs) were calculated for both 
the bulk and surface of the Si. One expects to observe 
similar DOSs for the bulk of Si and the deepest atom 
from the surface of the Si(lll). Our calculated surface 




(c) (d) 



FIG. 1: (Color online) (a) A Pb/Si(lll) slab symmetrically 
immersed in a vacuum to illustrate atomic arrangement of Pb 
adatoms (larger grey spheres) on 8 bilayers of Si(lll) atoms 
(smaller light spheres), (b) A schematic view of Pb overlayers 
to illustrate Tl (top), H3 (hep), T4 (fee) and bridge (inter- 
mediate point of the Tl and T4 or Tl and H3) sites of the 
Pb adatoms on the Si(lll) substrate. Top (c) and cross (d) 
views of the the most stable phase, Tl, including a bilayer of 
Si(lll) and a Pb overlayer located along the first layer of the 
Si(lll) bilayer. The (lxl)-unit cell is outlined in (b) and (c). 



DOS for Nsi < 8 predicts a wrong metallic behavior for 
the deepest Si atom. However, our result confirms that 
the DOS of the deepest Si layer approaches to the bulk 
DOS of the Si semiconductor by approaching Nsi to the 
8 bilayers. We avoid to present the DOSs of the substrate 
without the deposited Pb layers in this paper, as we just 
used them to ensure the validity of the Pb/Si(lll) thin 
films calculations. The Pb layers were then deposited 
over and below the prepared Si(lll) substrate to create 
a symmetric Pb/Si(lll) slab. In order to simulate the 
Pb/Si(lll) thin films, the created Pb/Si(lll) slab, as 
shown in Fig. [TJa), is symmetrically immersed in a vac- 
uum as well. The symmetric thin films causes to add 
inversion symmetry which gives rise to speed up the in- 
terface relaxation. The vacuum thickness is defined to 
be the distance between the top edge of the slab and the 
bottom edge of its next neighbor. In order to determine 
the vacuum thickness, total energy, work function and 
exerted forces on the surface atoms were calculated for 
the 8 bilayers Si(lll) structure versus various vacuum 
thicknesses. Our results, which are not presented here, 
show that 12 A is sufficient for the vacuum thickness to 
avoid interactions with the nearest neighbors of the slab 
along the Cartesian z axis. The total energy, work func- 
tion and forces were well converged in the vicinity of the 
obtained vacuum thickness - they did not show signif- 
icant changes by increasing the vacuum thickness more 
than the value obtained. All the Si underneath atoms 
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TABLE III: Relative total energies with respect to the en- 
ergy of the Tl phase, E(T1), in the presence and absence of 
the spin-orbit coupling (SOC). Here the E(T1) is arbitrarily 
chosen to be zero. 



SOC 


Tl (eV) 


T4 (eV) 


H3 (eV) 


No 


0.00 


0.41 


0.32 


Yes 


0.00 


0.25 


0.20 


No a 


0.00 


0.35 


0.18 



" -Reference [35 



and Pb adatoms are allowed to fully relax by adjusting 
their heights, but not their coordinates parallel to the 
Si(lll) surface. The top (Tl), fee (T4), hep (H3), and 
bridge (B2 - intermediate point of the Tl and T4 or Tl 
and H3) sites, as corresponding phases for the deposi- 
tion of the Pb layers on the Si(lll) surface, are shown 
in Fig. [TJb). In order to better visualize the most sta- 
ble phase, top and cross views of the (v^xv^-unit cell 
of the Tl phase are illustrated in Figs, [ljc) and (d), re- 
spectively. In this phase, as shown in Fig. [TJd), the Pb 
adatoms were positioned along the first layer of the the 
first Si(lll) bilayers. The (lxl)-unit cell as the inter- 
face structure for our full-potential all electron ab initio 
calculations is outlined in Figs.QJb) and (c). 



V. INTERFACE PROPERTIES 

A. Stability 

In this section, we would make sure about the most 
preferable configuration. Total energies were then calcu- 
lated for all of the illustrated structures in the Fig.QJb), 
apart from the bridge-site structure. The calculations 
were performed in the presence and absence of the spin- 
orbit coupling (SOC). Total energies relative to the Tl 
site are presented in Table. IIIII Our result with and 
without SOC shows that the Tl site is the most stable 
structure among the other configurations. This result is 
in agreement with experiment. [|| The energy, as shown in 
Table. IIIII increases from Tl site to T4 site and deceases 
from T4 site to H3 site. The later trend, as shown in 
Table. IIIII is not affected by the spin-orbit interactions. 
The result shows that the effect of spin-orbit interactions, 
keeping the trend, is to shift downward total energies of 
all the sites. Our results are also quantitatively com- 
parable with the pseudopotential calculations. [33j The 
energy of the bridge-site phase has not been calculated, 
since, as mentioned in the previous Sec. IIVI B2 is an in- 
termediate site of Tl and T4 sites or Tl and H3 sites. 
Therefore, one can most likely predict that the surface 
formation energy of the B2 site because of its symmetry 
is somewhere between Tl and T4 sites or between Tl and 
H3 sites. This prediction is in complete accord with the 



TABLE IV: Pb-Si bond length, d(Pb-Si), together with the 
method of calculations and exchange-correlation (XC) func- 
tional. 



Method 



XC-Potential 



d(Pb-Si)(A) 



FP-APW±lo 

FP-APW+lo 

FP-APW+lo 

PPW Q 

EXP. 6 



WC-GGA 
PBE-GGA 

LDA 

LDA 



2.72 
2.70 
2.68 
2.66 

2.66 ±0.03 < d < 2: 



±0.03 



" Reference |j. 
b Re ference\34 



pseudopotential results. [33| We will concentrate only on 
the lowest-energy Tl site from later on in the remaining 
subsequent sections. 



B. Si-Pb bond length 

Bond length of Pb-Si, d(Pb-Si), at the interface of 
our slab has been calculated for the most energetically 
stable phase, i.e., Tl site, as shown in Tab. IIIII The 
WC-GGA01, PBE-GGA[13 and LDA[26| functional 
were used for the exchange-correlation (XC) potential. 
Our full-potential results together with the pseudopo- 
tential results of the others [Jl as well as experimental 

The experimental data 



data[34j are listed in Tab. [TVj 
were obtained by surface X-ray diffraction based on the 
Pb/Si(lll)-(_? I) modelHI. This model contains 7 Pb 
adatoms in the (_f f)-unit cell [13]. Here we are only 
interested in the four of these lead atoms, since they are 
located in the top (Tl) site as our predicted most sta- 
ble phase (341]. The Pb-Si bond lengths for these 4 lead 
adatoms were measured (34l| to be between 2.66 ± 0.031 
and 2.98 ± 0.031.(13] The covalent radii of silicon and 
lead are R co (Si) = 1.181 and R co {Pb) = 1.471, re- 
spectively, and the bulk radius of lead is Rb u (Pb) = 
1.751. [34J Thereby the covalent bond lengths of Pb-Si 
have to be between R co (Si) + R co (Pb) = 2.651 and 
Rco{Si)+Rbu{Pb) = 2.93A, viz. 2.651 < d co (Pb-Si) < 
2. 93l [33] Our LDA full-potential result, as shown in 
Tab. IIVI is in excellent agreement with the LDA pseu- 
dopotential result of the others. [3] The LDA results, as 
shown in Tab. IIVI are very close to the strongest covalent 
bond length, i.e., 2.651. The result demonstrates that 
the PBE-GGA and WC-GGA give larger bond lengths 
which result in weaker Pb-Si bonds than LDA. This con- 
firms a trend shown in Tab. [I] in which the LDA un- 
derestimates but GGA overestimates lattice parameters 
compared to experiment. The weaker PBE- and WC- 
GGA bonds, as shown in Tab. IIVI are still in the covalent 
bond length interval, i.e. [2.65, 2.93 A]. The predicted 
covalent bonds by LDA and GGA show that the com- 
petition between the Pb-Si and Pb-Pb interactions gives 
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rise to saturate all the dangling bonds of the silicon sub- 
strate. Finally, it worths to mention that the trend seen 
in Tab. H] is reversed here, as shown in Tab. IIV( between 
PBE-GGA and WC-GGA in which PBE-GGA predicts 
lattice parameter larger than WC-GGA. In contrast to 
the lattice parameter trend, here, as shown in Tab. IIV[ 
the Pb-Si bond length within the WC-GGA is slightly 
larger than the Pb-Si bond length within the PBE-GGA. 
Normal and in-plane Pb-Si distances, viz. d z and d xy , 
were calculated employing pseudopotential method for 
the V7 x V3 mixed phase obtained from coadsorption 
of Pb and Sn on the Si(lll) surface. The Pb-Si bond 
lengths, d = (d xy + d 2 z ) * , were found to vary from 2.77 
A to 3.33 A depending on the positions of silicon atoms 
with respect to the lead adatoms for the (Pb, Sn)/Si(lll) 
V7 x phase. [l7t The Pb-Si bond length of the later 
\/l x V3 model with less in-plane Pb-Si distance can be 
more comparable with the bond length obtained from 
our top site phase. The less in-plane distance and its 
corresponding normal distance were calculated to be d xy 
= 1.79 A and d z = 2.11 A, respectively, which resulted 
in Pb-Si bond length d = (1.79 2 + 2.11 2 )^ = 2.77 A.[l7| 
The later bond length, 2.77 A, which is in the above men- 
tioned covalcnt bond length interval, i.e., [2.65, 2.93 A] 
and closer to our WC-GGA result, asserts that the Pb-Si 
forms a very strong bond.fllj 



C. Work function and surface formation energy 



TABLE V: The work function, cj> (eV), and surface forma- 
tion energy, E s (eV/A 2 ), per (lxl)-unit cell for the Pb at 
the interface (1 ML coverage) of the Pb/Si(lll) together with 
the method of calculations, kind of exchange-correlation (XC) 
functionals and interactions as well as existence or extinc- 
tion of the Si(lll) substrate. Since most of the presented 
pseudopotential results of the others in this table are for the 
hypothetical free-standing Pb(lll), our FP-APW+lo results 
without Si(lll) substrate are also included for the compari- 
son. 



Method 


XC-Pot. 


SOC Si(lll) 


<j> (eV) E s (eV/A 2 ) 


FP-APW+lo 


WC-GGA 


No 


Yes 


4.72 


-0.59 


FP-APW+lo 


WC-GGA 


Yes 


Yes 


4.65 


-0.61 


FP-APW+lo a 


WC-GGA 


No 


No 


4.23 


0.02 


PPW b 


LDA 


No 


Yes 




-0.47 


PPW C 


PW-GGA 


No 


No 


4.07 




PAW-PPW d 


GGA92 


No 


No 


3.83 


0.03 



"These free-standing Pb(lll) values are given for N=2. For 
N=l, the hypothetical free-standing film is so unstable that 
may no longer be considered as a physical system even 
within a pure theoretical study. 
b Reference % 
c Reference 3q 

d Reference\31l The given surface formation energy in the 
reference [33, 0.372 eV, is converted from eV to eV/A 2 . 
Hence the value of 0.372 eV is divided by the area of the 
(1 x l)-unit cell, 12.804 A 2 . 



We calculated the work function ,4>(eV), as the min- 
imum energy required to liberate an electron from the 
Fermi level (Ep) to a point with negligible kinetic energy 
at the center of the vacuum of the slab. The calculations 
were performed using the following formula: 



= E„ar. — Ef 



(1) 



where E vac is estimated by the averaged electrostatic 
Coulomb potential at the midpoint of the vacuum of the 
slab and Ep is the corresponding Fermi energy. The 
work function is a sensitive parameter to the surface con- 
ditions, since the liberated electron must move through 
the surface. The work functions were calculated at the in- 
terface of the slab with and without SOC and the results 
are presented in Tab. [Vj The result shows that the SOC 
causes to reduce the work function. The pseudopoten- 
tial results for the clean Pb(lll) surface [3a, [33] are also 
given in Tab. IVl within two different exchange-correlation 
functionals. The comparison shows that the various ver- 
sions of GGA can result in slightly different values. Since 
the pseudopotential work functions are available for the 
hypothetical free-standing Pb(lll), we also added our 
results without Si(lll) substrate in Tab. [Vl for the com- 
parison. Our free-standing work function is calculated for 
the second layer of the Pb(lll). A single lead layer with- 
out any substrate may not constitute a meaningful stable 
physical system. Our Pb(lll) work function which is in 



good agreement with the pseudopotential result would be 
compared with the results of our actual Pb/Si(lll) slab 
to elucidate the role of Si(lll) substrate. An overlook on 
the Tab. M shows that our calculated work functions tak- 
ing silicon substrate into account are in the same order 
of magnitudes when compared to the free-standing full- 
and/or pseudo-potential results. From the later point, it 
appears that the Si substrate does not significantly in- 
fluence the work functions. The later point confirms the 
observation of P. S. Kirchmann ct al.,[l3j], where theyJUJ 
found the agreement between their experimental QSE ob- 
servation for the actual Pb/Si(lll) and the pseudopoten- 
tial results 12j for the hypothetical freestanding Pb(lll). 

The surface formation energy per unit area at zero tem- 
perature is defined as:[35j 



Ef 



1 

2A 



(Eslab 



N S iE b s f k 



N Pb Ep^ k ) 



(2) 



where E s i a i, is the total energy of the slab, Nsi(Npb) 
and Egf (E h p^ k ) are the number and bulk energy of the 
Si (Pb) atoms in the unit cell, respectively. Here A is 
the area of the (lxl)-unit cell. The factor ^ is used 
in the above formula, since the slab has two surfaces, 
as shown in Fig. UK a), due to the embodied inversion 
symmetry. In principle one can use the above equation as 
the standard method to calculate the surface formation 
energy. However, it is well known that in practice the 
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surface formation energy employing Equ. [2l diverges. |38l 
H§| In order to overcome the divergence problem of the 
standard Equ. [2 in this paper we have used the following 
formula[38j to calculate the surface formation energy: 



E f = —{Esiab-Nsi&E, 



(3) 



where AEg 



Eg 



Eg- 1 and A£& = E*L - E?' 1 



Pb ■ 

We have already used Equ. [3j with Npf, = 0, in the pre- 
vious Sec. IIVI to calculate the surface formation energy 
of the pure Si(lll) substrate. Our calculated surface 
formation energies in the presence and absence of the 
spin-orbit coupling (SOC) along with the pseudopoten- 
tial results for both of the actual Pb/Si(lll) interface^ 
and hypothetical clean Pb(lll) surface [37| are given in 
Tab. |Vl Our result shows that the effect of SOC is small 
as well. Spin-orbit (SO) coupling breaks down the space 
inversion symmetry by imposing a preferable direction, 
which gives rise to spin splitting. The inversion symme- 
try can be also broken at an interface or a surface. The 
Rashba-effect[!o| on metallic surfaces is a phenomenon, 
which originates both from spin-orbit coupling and the 
lack of inversion symmetry at surfaces. [41| A Rashba- 
type spin-orbit splitting with no coverage dependency 
was found for quantum well states formed in ultrathin 
Pb films on Si (111). [3 Hugo Dil et al.0 found for the 
states at 0.15 and 0.4 eV this splitting to be 14 and 15 
meV, respectively. We found the same splitting in the 
vicinity of Fermi level for the valence bands. The SO 
coupling within this small splitting, as shown in Tab. [Vj 
causes to reduce the work function and surface forma- 
tion energy by 0.07 eV and 0.02 eV/A 2 , respectively. As 
shown in Tab. [Vj the surface formation energy is a neg- 
ative value taking substrate into account, i.e., for the 
actual Pb/Si(lll) case. On the contrary it is a tiny pos- 
itive value without substrate, i.e., for the hypothetical 
freestanding Pb(lll) case. The obtained negative value 
when compared with the tiny positive value confirms that 
the Pb/Si(lll) is energetically more favorable to form 
that particular structure than the clean Pb(lll) surface. 
The later point is in complete accord with the pseudopo- 
tcntial results, as can be seen from Tab. [V] Therefore, 
the effect of substrate is obviously to make the Pb layers 
energetically stable. 



D. Density of states (DOS) 

Total and partial densities of states (DOSs) were cal- 
culated at the interface for both of the underneath Si 
and the ovcrlayer Pb atoms. The calculations were per- 
formed in the absence and presence of the spin-orbit cou- 
pling (SOC) for the slab shown in Fig. [Ha). Total DOSs 
with and without SOC, which are not shown here, were 
compared. Despite the large atomic number of Pb, the 
comparison shows that the effect of SOC is not very sig- 
nificant on the valence electronic structure and can be 
ignored. This is in accordance with the experiment in 
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FIG. 2: (Color online) Partial densities of states for (a) s Pb, 
(b) p Pb, (c) p x+y Pb, (d) p 2 Pb, (e) s Si, (f) p Si, (g) p x+y 
Si, (h) Pz Si. 



which the magnitude of the Rashba-type spin-orbit split- 
ting were found for this thin films to be so small that 
can be detected by the spin integrated angle resolved 
photoemission (ARPES) method. [16| The spin splitting 
originates from all of the layers and thereby should be 
calculated by integrating over the entire layers. [l|| The 
spin splitting is small in Pb/Si(lll) ultrathin films, be- 
cause the signs of the splittings of the Pb-vacuum and 
Si(lll)-Pb interfaces are opposite in direction due to the 
larger phase shift at the Pb to the vacuum interface than 
the phase shift at the Pb to Si(lll) interface. [l(| Thus 
the total splitting which partly cancels each other results 
in an ignorable spin-orbit coupling effect in Pb/Si(lll). 
Therefore from now on we do not include the spin-orbit 
interactions in all of our surface calculations. Partial s 
and p DOSs are shown in Figs.^a) and (b) for Pb and in 
Figs.[2]Jc) and (f) for Si atoms at the interface of the slab, 
respectively. The result shows that the s states of the Pb 
and Si atoms are almost distributed over the energy in- 
terval of (-10, -5) eV, while their p states are nearly in 
the (-4, Ef = 0) eV energy window. This demonstrates 
that the s states are lower in energy than the p states 
for both of these Pb and Si atoms. Consequently, the s 
states of Pb and Si are more localized than their p states, 
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which are closer to the Fermi level when compared with 
the former s sates. Since the p orbital manifest itself 
more closely to the Fermi level, the p sates show more 
itinerant character than the s sates. Hence, the p elec- 
trons play more important role than the s electrons in 
the Pb-Si bonding. Now we can discuss the kinds and 
qualitatively the strengths of the Pb-Si bonds, keeping 
three points in mind: (i) s (p) electrons are more local- 
ized (itinerant) than p (s) electrons for both of the Pb 
and the Si atoms, and (ii) below the Fermi level inter- 
val energy (-10, -5) eV of s states do not overlap interval 
energy (-4, 0) eV of p states for both of the Pb and Si 
atoms, as well as (iii) Pb sits on-top of Si site vertically 
along z direction, as shown in the top site configuration 
in Fig. [1] (b) . The result shows that the 6s states of Pb 
adatoms are weakly hybridized with the 3s states of the 
Si substrate. According to the above point (ii) the en- 
ergy window of 3s Si well overlaps the energy window of 
6s Pb, as shown in Figs. [DJa) and (e). One may then 
promptly speculate that the 6s Pb states can be well hy- 
bridized with the 3s Si states. However, according to the 
point (i) both of the 6s and 3s states are nearly local- 
ized and so do not substantially contribute in bonding. 
Thereby the overlap of the 3s Si and 6p Pb DOSs can 
only result in an extremely weak er-likc bond, which may 
be due to its weakness neglected. The 6s Pb states are 
even much more weakly bonded to the p Si states com- 
pared to the former extremely weak cr-like Pb-Si bond. 
According to the property (ii) energy windows of 6s Pb 
and 3s Si do not overlap, because the s-Pb states are 
over (-10, -5 eV), while the p-Si state are in (-4, eV), 
as can be clearly seen from Figs. [2ja) and (f). This in 
turn gives rise to even weaker bond between 6s Pb and 
p Si states than the last 6s Pb and 3s Si bond. This is 
also the case, as shown in Figs, [^b) and (e), for the p 
Pb and s Si DOSs. Till now, there are no any evidence 
of the appearance of Pb-Si bonding. Let us then discuss 
the bonding states by concentrating solely on the most 
important p states of the Pb and Si. First, the result 
shows that the energy window of p Pb well overlaps the 
energy window of p Si, as shown in Figs. [UJb) and (f). 
Second, according to the above mentioned (i) point, the 
p Pb and Si electrons are more itinerant than their s elec- 
trons. Third, there are a lot of p Pb and p Si states at 
an energy very close to the Fermi level, see two sharp 
peaks in Figs. [5£b) and (f). Therefore, the p states of 
the Pb and Si atoms can be well hybridized with each 
other. In order to determine the kind of Pb-Si bond, the 
p states of Pb and Si were further decomposed, as shown 
in the right panel of the Figs. [H to their partial p z and 
p xy states. The distributions of the p z and p xy DOSs of 
both Pb and Si atoms show that they are completely hy- 
bridized with each other. Here, the z axis is taken to be 
perpendicular on the surface and x,y are parallel to the 
surface. The Pb overlayer atom is just located above the 
Si underneath atom along the z-Cartesian coordinate in 
the top site configuration, as shown in Figs.[ljb). Hence, 
one lobe of the p z orbital of Si underneath atom over- 



laps with one lobe of the p z orbital of Pb overlayer atom. 
Consequently, hybridization of p z orbitals of Pb and Si 
results in a cr-bond. The strength of the p z Pb and p z 
Si cr-bond can be estimated by observing that the two 
sharp peaks at an almost single energy nearby the Fermi 
level originate from the p z Pb and p z Si DOSs, as shown 
in Figs. [IJd) and (h). This can be taken as an indica- 
tion to the fact that there is a tightly cr-bond between 
the the p z Pb and p z Si at the interface. This result is 
in excellent agreement with the pseudopotential result of 
the others[I3 for Pb/Si(lll) V3 x y/7 phase, in which a 
strong hybridization between the P z Pb orbital and the 
P z dangeling bonds of the substrate had been reported. 
As shown in Figs. &[c) and (g), two lobes of the p xy Pb 
overlap with two lobes of the p xy Si as well, which gives 
rise to a 7r-bond. The strength of the 7r-bond cannot be 
so large, because the p xy Pb and p xy Si DOSs are broad- 
ened, as shown in Figs.[Hd) and (h), in contrast to the 
sharp narrow p z Pb and Si DOSs. Thereby, a feeble ir- 
bond is constituted by the the p xy Pb and p xy Si at the 
interface. We conclude this discussion by stating the fol- 
lowing points. There is not a pure n or a state in the 
Pb-Si bond. Instead there is a mixed state composed of 
7T and cr states in the Pb-Si bond. The combination of 
these two feeble tt and strong cr bonds can result in a 
strong covalent Pb-Si bond. The later point is consistent 
with our last discussion in Sec. IV Bl regarding the bond 
length of Pb-Si. In summary, the Pb overlayer can be 
strongly absorbed by the Si substrate. 



VI. THIN FILMS PROPERTIES 

A. Total Energy and Energy Differences 

Total energy per number of layers, E(N)/N, is shown 
in Fig. [3] (a) as a function of the film thickness, N. The 
result shows that the E(N) /N increases by N and asymp- 
totically approaches to a constant value. The asymptotic 
behavior can be realized, if we notice that the motion 
of Pb electrons is confined in one direction while it re- 
mains free to move on the surface. The later point itself 
is due to the fact that the Pb/Si(lll) is a prototype of 
the metal-on-semiconductor system with a high Schottky 
barrier againstpenetrating valence electrons of Pb into 
the substrate. [43. |43| The confinement in the nanoscale 
causes wave vector k to quantize in a perpendicular di- 
rection on the surface, which results in discrete energy 
levels associated with the so-called quantum-well (QW) 
states. The simplest way to describe this confinement is 
then by the one dimensional QW.[44] In a direction nor- 
mal on the surface there is no translation symmetry. The 
potential in the xy plane is periodic and potential in the 
z direction represents the QW. The energy levels in the 
free electron approximation are given by: 
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FIG. 3: (Color online) (a)Total energy per layer (solid line) 
and energy differences (dotted line) of Pb/Si(lll) versus num- 
ber of monolayer (ML) coverages, N. (b)Surface formation 
energy (dotted line) per (lxl)-unit cell and work Function 
(solid line) of Pb/Si(lll) versus number of monolayer (ML) 



coverages, N. 



where m is the free electron mass and d is the thickness 
of the film. The above equation clearly shows that the 
total energy approaches to a constant value by increasing 
the film thickness. This equation elucidates recognition 
of the asymptotic behavior shown in Fig. [3] (a). Mono- 



layer per atom energies 



E(N) 
N 



were calculated for the free- 



standing Pb(lll) slab as a function of monolayers N.fToj 
In agreement with our above discussed asymptotic be- 
havior for the actual case of Pb/Si(lll), it was shown 
that »r of the hypothetical freestanding Pb(lll) sys- 



tem gradually approaches a constant value. [Io| 

The energy difference between two successive layers is 
defined as: 



5E(N) = E(N)-E(N - 1). 



(5) 



The SE(N) versus number of layers, N, is also shown in 
Fig. EJa). The result shows that the energy difference 
oscillates as a function of the number of monolayer (ML) 
coverages with the period of A = 2ML. The oscillation 
is attributed to the quantum size effect (QSE). Similar 
QSEs were detected by the pseudopotcntial calculations 
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FIG. 4: (Color online) Surface formation energy (dotted line) 
per (1 x l)-unit cell and work Function (solid line) of freestand- 
ing Pb(lll) versus number of monolayer (ML) coverages, N. 



for the freestanding Pb(lll) thin films. [Icl] One expects 
that the underlying physics behind of the observed QSE 
in the 6E(N) quantity, as shown in Fig. [3Ja), can be 
related to the electronic structure of the system. We 
would postpone this point to a more appropriate time 
in Sec. IVI C[ where the electric field gradient (EFG) as 
an extremely sensitive quantity to the valence electron 
charge density distribution is discussed. 



B. Work function and surface formation energy 

The work function, as discussed in the previous 
Sec. IV C[ has been calculated utilizing Equ. [T] The re- 
sult is presented in Fig. Efb) as a function of Pb thin 
films thickness for the actual Pb/Si(lll) slab. The result 
shows an oscillatory behavior of the work function ver- 
sus the film thickness with the period of A = 2ML. The 
oscillation is quickly damped and the work function is go- 
ing to be converged through a few number of Pb layers. 
The work function oscillation was previously attributed 
to the quantum size effects (QSE) and predicted[l2j by 
using a pseudopotential calculation for the case of a free- 
standing Pb slab. The pseudopotential result for the 
hypothetical clean Pb(lll) thin films shows that the os- 
cillation is much more slowly damped and the work func- 
tion is converged at a very large number of Pb layers 
when compared with our full-potential results for the 
real Pb/Si(lll) case. The different oscillatory behav- 
ior might be related to the effect of the Si substrate. But 
one has to notice that the discrepancy may come from 
the two different approaches. Thus we calculated the 
work function for the clean Pb(lll) surface within our 
full-potential method as well. The result as shown in 
Fig. 2] shows that our full-potential work functions are 
also no longer rapidly damped for the clean Pb(lll) case 



9 



in agreement with the pseudopotcntial result. This ver- 
ifies that the oscillatory behavior discrepancy indeed re- 
sults mainly from the existence and extinction of the Si 
substrate and not from the method of calculations. This 
authenticates that the occurred QSE in the work func- 
tion is influenced by that of Si(lll). The effect of Si 
substrate on the work function confirms the calculation 
of Dil et al.,[l4[ where thev(l4j similarly found that unlike 
for Pb on graphite, the Pb overlayer lattice structure is 
influenced by the Si(lll). The effect of our Si(lll) sub- 
strate on both of the geometry and electronic structures 
of the Pb overlayer atoms is also in agreement with the 
influence of the Cu(lll) substrate on the geometry struc- 
ture of the Pb(lll) layers reported by Materzanini and 
coworkers. [Tpj In the later work, the effect of substrate 
was studied in a tricky manner without including directly 
Cu(lll) underneath atoms by compressing the super- 
cell dimensions in the surface plane by 3.3%. [l(| They 
found that the explicit inclusion of the Cu(lll) substrate 
would be highly desirable to explain the remaining quan- 
titative differences between theory and experiment . [Toj 
It seems that the later result contradicts our last find- 
ing in Sec. IV CI We found in Sec. IV CI that the order 
of magnitude of the work function could not be signif- 
icantly affected by the substrate in agreement with the 
observation presented in Ref. [l3|. The contradiction can 
be resolved if we notice that the work function is not so 
sensitive that its order of magnetite, but its oscillatory 
behavior, can be affected by the substrate. The electric 
field gradient (EFG) is much more sensitive quantity than 
the work function to the surface states. We shall show in 
subsequent Sec. IVI Cl that the EFG is so sensitive to the 
substrate that not only its oscillatory behavior but also 
its value can be affected by the Si substrate. 

We calculated the surface formation energy versus the 
film thickness, N, based on the Equ. [2] The result, which 
is not shown here, varies linearly with respect to N. The 
linear variation is consistent with the first scheme de- 
scribed in Ref. [H. We then recalculated the surface for- 
mation energy versus N according to the Equ. [3] The 
result, as shown in Fig. EJb), manifests a deviation from 
the linear behavior on going from fourth to the fifth layer. 
The deviation is consistent with the third scheme dis- 
cussed in Ref. [3a. Here we did not use the fourth scheme 
given in Ref. [38|. According to the fourth scheme the Bulk 
energies in the standard Equ. [2] can be found by taking 
the slope of a fitted straight line to all of the slab to- 
tal energies versus N. It is generally believedJH, [H, Hfl 
that the divergence problem disappears employing the 
fourth scheme for smaller N. However, in this paper we 
aim to elucidate the role of substrate. The goal of this 
paper can be achieved even by the third scheme which is 
identical to the the Equ. [31 Therefore we calculated the 
surface formation energy versus N for the clean Pb(lll) 
surface. The result is shown in Fig [4] which can be 
compared with the surface formation energy of the ac- 
tual Pb/Si(lll) slab as shown in Fig. |3Jb). Furthermore 
the clean Pb(lll) surface formation energy behaves com- 
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FIG. 5: (Color online) The main component of the electric 
field gradient, V zz , versus number of monolayer (ML) cover- 
ages, N. 



pletely different from that of its supported slab by the Si 
substrate. The different behavior demonstrates the effect 
of Si on the surface formation energy. 



C. Electric Field Gradient (EFG) 

The electric field gradients (EFGs) [H, H3] were calcu- 
lated at the Pb sites for our slab. The main component 
of the EFG tensor, V ZZl is presented as a function of film 
thickness, N, in Fig. [5] The sign of the V zz is negative, as 
shown in Fig. O for all the layers. The result shows that 
for 1 monolayer coverage, i.e., N = 1, the absolute value 
of the V zz is a very large number. As shown in Fig.0 the 
| V zz | is drastically reduced on going from first to the sec- 
ond layer. Despite the reduction is significant, the \V ZZ \ 
values are still considerably large for N > 2. Our result 
shows that for the larger N the V zz oscillates with re- 
spect to the number of monolayer (ML) coverages. The 
period of oscillation here is also A = 2ML. Therefore 
our ab initio calculations clearly elucidate the quantum 
size effect (QSE) in the electric field gradient (EFG). In 
order to interpret the above observation shown in Fig. [5] 
a cartoonlike electron charge distribution is illustrated in 
Figs.[H The electric field gradient (EFG) is an extremely 
sensitive quantity to the deviation from t he sp herical va- 
lence electron charge density distribution. [47l. lil. |49|. |50| 
The EFG, as shown in Fig. [6] (b), is zero for a spherical 
charge density distribution, which is the case for the cu- 
bic or higher point group symmetry. The EFG, as shown 
in Figs. [6] (a) and (c), can be nonzero at an atomic site 
depending on its point group symmetry for an asphcri- 
cal charge density distribution, which is the case for our 
Pb sites. The EFG can be calculated within the APW 
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FIG. 6: (Color online) A schematic representation of elec- 
tron charge density distribution (ECDD). (a) A longitudinal 
ECDD deformation can yield a nonzero negative value for 
the main component of the electric field gradient, V zz , for 
angular 1=2 and azimuthal m=0 quantum numbers, (b) An 
undeformed spherical ECDD yields a zero value for the V zz , 
for every I and m. (c) A transverse ECDD deformation can 
yield a nonzero positive value for the V zz , with 1=2 and m=0. 



method by V zz oc f[p(r)Y w ]/r 3 dr.^ The EFG can be 
nonzero for an aspherical charge density distribution and 
obtained directly from the I = 2, m = component of 
the potential expansion inside the muffin-tin spheres. [50j 
The sign of the EFG depends on the orientation of the 
anisotropic charge density distribution. The sign is de- 
fined to be negative (positive) for the charge distribution 
shown in Fig. [5] (a) (Fig. [5] (c)). Our surface Pb atoms 
for 1 monolayer (1 ML) coverage are only supported by 
their underneath Si(lll) substrate atoms, since for N = 
1 there is no more Pb overlayer on the first deposited Pb 
layer. The Pb-Pb bonds in the (x, y)-plane compensate 
each other. However, there is no Pb-Pb bonds along the 
z Cartesian direction to compete with the Pb-Si bonds 
to compensate them. The Pb-Si bond, as discussed in 
Sec. IV Bl is a strong covalent bond. The uncompensated 
covalent Pb-Si bond results in a dramatic electron charge 
density deformation along the z direction. Consequently 
one can, according to the Fig. [5] (a), first expect that 
the sign of the EFG is negative and second its value is 
dramatically large. These expectations are in excellent 
agreement with our result as shown in Fig. [5] For 2 ML 
coverages the second Pb layer is supported by the first 
Pb layer. The Pb-Pb bonds, which are a gain compen- 
sated in the (x, y)-plane, are much weaker 3J| than the 
Pb-Si bonds along the z Cartesian direction. The weaker 
bonds results in a less charge density deformation which 
gives rise to a less value of the EFG at the Pb site for the 



second layer. The later point clearly explains the signifi- 
cant drop of the EFG on going from N = 1 to 2 as shown 
in Fig. Although the EFG is drastically reduced by 
adding the second layer, it can be still large due to the 
anisotropy of the Pb surface. The anisotropy comes from 
the fact that for N = 2 there is no Pb-Pb bond over the 
second layer to compensate the Pb-Pb bonds between the 
first and second layers. 

Here we would back to the postponed point to clarify 
that the source of the quantum size effect (QSE) origi- 
nates from the electronic structure of the system. In or- 
der to reach the purpose one can first make a connection 
between the EFG and the density of states (DOS). Fortu- 
nately an astonishing connection was previously made [52j 
between the EFG and the total DOS at the E F , N F . 
There it was found an approximately linear relation be- 
tween the main component of the electric field gradient 
and the total density of states at the Fermi level, viz. 
Vzz « Nf, in the Celn^ compound. [52j Therefore from 
one side there can be a relation between the V zz and the 
Np in our thin films. On the other hand here we have 
shown that the V zz oscillates as a function of film thick- 
ness. Consequently, there might be a relation between 
QSE and the electronic structure, or more specifically, 
between QSE and the total DOS at the Fermi level (Np). 
Our finding regarding the later relation between QSE and 
Np verifies the previously reported relation between Np 
and the electron spillage length into the vacuum and as 
a result between Np and the step height of the layer 
This illuminates the underlying physics behind of 
the quantum size effects (QSEs). Our EFG study shows 
that the oscillations of a sensitive physical quantity with 
respect to the film thickness originates from the beat- 
ing of the valence electron charge density distributions 
(VECDD). The VECDD deformation can change layer 
by layer depending on the bonds strength. The effect 
of the Si substrate, which is shown to be of significant 
importance in the strength of the bonds, plays an impor- 
tant role in the QSE. One suspects that the freestanding 
model can be applicable for those kinds of the thin films 
that the strength of their substrate-overlayer bonds are 
comparable with the strength of their overlayer-overlayer 
bonds. The later point is not the case for our thin films, 
since the strengths of the Pb-Si and Pb-Pb bonds are not 
comparable. [34| In conclusion our QSE calculation shows 
that our system is not an ideal paradigm to freestanding 
films. 



VII. CONCLUSIONS 

The effects of Si(lll) substrate were investigated on 
the physical properties of the Pb/Si(lll) thin films. The 
investigations were performed within the density func- 
tional theory (DFT) employing the augmented plane 
waves plus local orbital method (APW+lo). We used the 
PBE-GGA and WC-GGA for the exchange-correlation 
functional. Our result shows that the WC-GGA is more 
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reliable to study this system. We included spin-orbit in- 
teractions in our the Kohn-Sham Hamiltonian. Our re- 
sult shows that the effect of spin-orbit coupling is not of 
significant importance for this thin film. Several struc- 
tures were considered for the slab. Our ab initio full- 
potential calculation shows that the top site (Tl) is the 
most stable phase. The electronic structures at the inter- 
face of the most stable Tl phase were studied. The study 
shows that the Pb and Si at the interface are strongly 
bonded by a and weakly by tt bonds. The strong cova- 
lent Pb-Si bond is consistent with experiment. In order to 
elucidate the role of Si substrate the Quantum size effects 
(QSEs) were studied for the Tl configuration. The study 
is carried out by calculating the work function and sur- 
face formation energy as well as the electric field gradient 
(EFG) at the Pb sites. The calculations were performed 
as functions of the thin film thickness (N). We connected 
the oscillatory behavior of our calculated physical quan- 
tities to the electronic structures at various layers of the 
system. The connection has been made by the concept of 
the electric field gradient (EFG) as an extremely sensitive 



quantity to the electronic structure. The EFG provides 
a reliable approach to undertake the underlying physics 
behind of the QSE. Therefore the electric hypcrfinc inter- 
action has merit as a measurement method to be further 
used in thin films. Our result shows that the effect of 
Si substrate depends on the sensitivity of the physical 
quantities to the valence electron charge densities. The 
Si substrate can be considerable for the sensitive quan- 
tities such as electric filed gradients (EFG), while it can 
be ignored for less sensitive quantities such as the value 
of the work function. 
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